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Abstract —In a mobile ad hoc network (MANET), effective 
prediction of time-varying interferences can enable adaptive 
transmission designs and therefore improve the communication 
performance. This paper investigates interference prediction in 
MANETs with a finite number of nodes by proposing and using 
a general-order linear model for node mobility. The proposed 
mobility model can well approximate node dynamics of practical 
MANETs. In contrast to previous studies on interference statis¬ 
tics, we are able through this model to give a best estimate of 
the time-varying interference at any time rather than long-term 
average effects. Specifically, we propose a compound Gaussian 
point process functional as a general framework to obtain 
analytical results on the mean value and moment-generating 
function of the interference prediction. With a series form of 
this functional, we give the necessary and sufficient condition 
for when the prediction is essentially equivalent to that from a 
Binomial Point Process (BPP) network in the limit as time goes 
to infinity. These conditions permit one to rigorously determine 
when the commonly used BPP approximations are valid. Finally, 
our simulation results corroborate the effectiveness and accuracy 
of the analytical results on interference prediction and also show 
the advantages of our method in dealing with complex mobilities. 

Index Terms —Interference Prediction, General-order Mobility, 
Compound Gaussian Point Process Functional, Gaussian BPP, 
Mobile Ad hoc Networks. 

I. Introduction 
A. Motivation and Related Work 

I N MANETs interference plays a pivotal role, through the 
Signal to Interference Ratio (SIR), in contributing to the 
Quality of Service (QoS). In contrast to static wireless net¬ 
works, the distances between interferers and receiver change 
dynamically during the times of communications because of 
the mobilities of interferers and receiver. As a result, the re¬ 
ceived signal is affected by fluctuating interferences generated 
by these mobile nodes. 

Interference analysis can help to discover and exploit the 
regularity of these time-vary interferences. For example, it 
helps to understand the statistical performance of communi¬ 
cation under such interference, e.g., outage probability. Due 
to disturbances on nodes’ movement or our incomplete infor¬ 
mation of the node locations, the trajectories of these nodes 
are often accompanied with uncertainties when analyzing 
interference in MANETs. Stochastic Geometry 0- 0 is a 
powerful tool for describing the random pattern of mobile 
nodes as a Point Process (PP) at each time instant. For 
mobile networks with high mobility nodes, the authors in |3j 
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analyzed the local delay, which is a functional of the inter¬ 
ference Moment Generating Function (MGF), by modeling 
node locations as a Poisson Point Process (PPP). Indeed, if 
one just focuses on one time instant, the existing results on 
interference analysis in static networks can be directly used 
in highly mobile networks (e.g., for PPP networks 0-0 or 
for Binomial Point Process (BPP) networks Q, @). Highly 
mobile networks imply that the node locations at two different 
time instants have almost no correlation, if the node velocities 
are sufficiently large. Nevertheless, for most practical cases, 
we are more interested in interference analysis in MANETs 
with finite nodal velocities. 

In order to analyze interference in MANETs with finite 
node velocities, mobility models are required to capture the 
node locations at each time instant. A summary of mobility 
models and their corresponding point processes in prior studies 
on interference statistics of MANETs is provided in Table [I] 
By employing random walk. Brownian motion and random 
waypoint model^J the statistics of interference were analyzed 
in m-m- For random walk and Brownian motion models, 
the approximate distribution of aggregate interference was 
given in ED- p~8) . The mean of the aggregate interference 
was analyzed in |17), and upper bounds in time-correlation 
for aggregate interference and outage probability were given 
in (l6j and ijTTJ . For the random waypoint model, similarly, the 
approximated distribution of aggregate interference was given 
in 03- 0D> and the mean of the aggregate interference was 
analyzed in [ [17] . 

However, the existing analysis methods only considered 
special or limiting scenarios where statistics of interference 
are identical at every time instant. For the random walk 
and Brownian motion models in 03 - GD- the means or 
pdfs (probability density functions) of aggregate interference 
at every time instant are the same due to the assumption 
that the initial node distribution is uniform. Consequently, 
node dynamics do not provide any useful information in the 
interference analysis. For the random waypoint model in [ 171 ], 
09’ the idea was to wait an infinite time for the node 
distribution to converge to one limiting distribution that can be 
regarded as a PPP with a quadratic polynomial intensity OS- 
Based on that kind of PPP, mean, outage and approximated 
distribution for aggregate interference were given. As a result, 
not only do node dynamics not contribute to the analysis, but 
also the results are only valid in infinite time. It should be 
noted that in many applications what we would like to know 

1 The constrained i.i.d. mobility model in 1151 is similar to the highly mobile 
network, thus we mainly discuss the interference analysis in MANETs under 
random walk, Brownian motion and random waypoint. 


TABLE I 

Summary of Mobility Models and Their Corresponding Point Processes in Prior Studies on Interference 

Statistics of MANETs 
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Reference 



Mobility Model 


Point Process 


Highly mobility networks 
Highly mobility networks with 
randomly actived interferes 
Highly mobility networks 
Constrained i.i.d. mobility 
Random walk 
Brownian motion 
Random waypoint 


PPP (time independent) 

PPP (time homogenous for sufficient large time) 

BPP (time independent) 

PPP (time homogenous, non-Markov 117|) 

PPP (time homogenous) 

PPP (time homogenous) 

PPP (time homogenous) with quadratic polynomial intensity 


Time independent: the locations of nodes change independently from one time instant to the next. 

Time homogenous: the locations of nodes are correlated between different time instants but have the same pdf. 


are the statistics of aggregate interference within a finite time 
window rather than an infinite one. 

Furthermore, the existing analysis on interference between 
multiple time instants typically focused on the time corre¬ 
lation (9) 11 11. p6) , (17j, which provides no further infor¬ 
mation beyond a linear relation. To design a communication 
strategy (e.g., transmission power control), we are interested 
in knowing and exploiting the statistics of interference at a 
future time of interest, which cannot be derived from a simple 
time correlation. We assert that interference prediction can 
provide us more effective information than time correlation. 
Beyond the simple temporal correlation, fl4| proposed a joint 
temporal characteristic function of interference for multiple 
time instants, and 0, (12), (B) investigated the condi¬ 
tional/joint outage/success probabilities over time. Despite the 
comprehensive characterizations of the temporal statistics, the 
studies in d, d assumed networks with either fixed 
or independent node locations from one time slot to the next. 
Hence, their results cannot be used for interference prediction 
in realistic network with practical mobility models. 

In fact, another limitation of current studies on interference 
in MANETs is the mobility model. For example, the random 
walk. Brownian motion and random waypoint models are often 
inadequate to describe many kinds of mobilities in the real 
world, like mobilities constrained by the physical laws of 
acceleration and velocity (20). Modelling should include all 
kinds of communicating objects that have the ability to move. 
With the development of automation, the need for communi¬ 
cation between robots is increasing. For civil use, unmanned 
aircrafts offer new ways for commercial enterprises and public 
operators to increase operational efficiency, decrease costs, and 
enhance safety ED- The node mobilities in such scenarios 
have complex dynamics and cannot be captured by random 
walk. Brownian motion, random waypoint mobility models 
or even the Gauss-Markov mobility model (20) , (22) that 
takes acceleration and velocity into account. Therefore, it 
is desirable to develop a more general mobility model that 
is capable to describe mobile nodes governed by complex 
mobility dynamics, and then use it for interference prediction. 

B. Our Contributions 

In this work, we focus on interference prediction in 
MANETs with general mobilities having a finite number of 


nodes. Compared to time correlation, interference prediction 
can give more effective information, i.e., providing the best 
estimate of the interference level at a future time instant based 
on the knowledge at the current time. By developing a general 
mobility model, the predictions can be used in a wide range 
of MANETs. 

The main contributions of this work are: 

• We propose a General-order Linear Continuous-time 
(GLC) mobility model to describe the dynamics of mov¬ 
ing nodes in practical applications. The random walk. 
Brownian motion and Gauss-Markov mobility models 
in 0, fig, (20) can be regarded as special cases of 
the GLC mobility model with discretizations. In this 
framework, the random walk and Brownian motion turn 
out to be first-order linear mobility models, and Gauss- 
Markov model is a second-order linear mobility model. 

• Based on the GLC mobility model, the mean and 
Moment-Generating Function (MGF) of interference pre¬ 
diction on a mobile reference point at any finite time into 
the future are derived and analyzed. In order to simplify 
the expression for mean and MGF under different path 
loss functions and multipath fading, a Compound Gaus¬ 
sian Point Process Functional (CGPPF) is defined and 
expressed in a series form. 

• Apart from interference prediction at finite time into 
the future, we also give the necessary and sufficient 
condition for when our predictions converge to those 
from a Gaussian BPP as time goes to infinity. This result 
provides a guideline on when the previous studies on 
interference statistics with BPP are relevant. 


C. Paper Organization 

In Section [II] we present the case of Uniform Circular 
Motion (UCM) as an example to motivate our general model. 
In Section [HI] we define and analyze the GLC mobility model. 
Additionally, the dynamic reference point is defined to study 
the relative node locations of interferers relative to the mobile 
receiver. In Section [IV] the mean and MGF of the interference 
prediction are analyzed. The numerical examples are given in 
Section [V] to illustrate the effectiveness of our approach and 
corroborate our analytical results. 
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(b) 


the mobility model can be approximated by UCM in the 2-D 
plane as two coupled differential equations 


= u>ix\ 2) (t) +w ( i 1 \t) 
xf\t) = + wf\t) 


where the subscript i denotes the /th node, and i = 
1,2 ,N, where N is the number of nodes in the MANET. 


{t)) is the location of node i in the 2D plane at 
time t, and w* is the UCM angular velocity, and the initial 
vector Xj(fo) = [x^ito), x[ 2 \to)] T determines the initial 
location and radius of node i. w^\t) and w^ 2 \t) stand for 
additive disturbance due to airflow. 

Consider a duration of time in which UAV 1 and 2 com¬ 
municate with UAV j and k, respectively. For UAV 0, it tries 
to receive information from UAV q. Therefore, there are two 
interferes that affect the signal reception at UAV 0. 


„( 2 ) 


The aggregate interference at UAV 0 is shown in Fig. 1(b) 


from which we see that the received aggregate interference at 
UAV 0 is periodically changing and does not converge to a 
constant value independent of time t. Interference predictions 
should be adaptive to the node dynamics and therefore make 
use of the characteristics of mobility models. 


III. General-order Linear Mobility Model 

In this section, the General-order Linear Continuous-time 
(GLC) Mobility Model is proposed and the corresponding 
statistics of node distribution is given. 


Fig. 1. (a) Three UAVs scan a target by using vision sensors and share 

information through a wireless network, (b) The received interference UAV 
0 versus time t. 


A. General-order Linear Mobility Model and Node Distribu 
tion 


D. Notation 


E[ ] denotes the mean of a random variable or a random 
vector, and E a [-] denotes the expectation operator with respect 
to a. Cov[ ] denotes the covariance of a random vector, and 
x(t) = dx(t)/dt denotes the derivative of stochastic process 
x(t) at time t in the mean square sense 1231. 


Consider a network having N nodes in d-dimensional space 
(e.g., d = 2 means nodes move in a 2D plane). For each node 
i, where i £ {1, 2,..., TV}, we model it employing the state- 
space model with additive uncertainties given by the stochastic 
differential and algebraic equations [24| 

ij (t) = A iXi (t) + w i ( t ) (2) 

y i(t) = CjXj(f), (3) 


II. Interferences in MANETs with Uniform 
Circular Motions 

In this section, we investigate the interference that arises in 
nodes exhibiting Uniform Circular Motion (UCM) to model 
UAVs carrying out a scanning task. This simple model will be 
seen to be a special case of the general GLC mobility model 
that we develop in Section m 

Consider N £ Z + UAVs carrying out a scanning task for 
a target, see Fig. |l(a)| which has many applications such 
as rescue operations, monitoring applications, etc. This target 
(the central pink circle) could be a point or an area. Due to 
limitations of their visual fields, each UAV can only acquire 
partial information about the target. Therefore, in order to 
have a better understanding of the target, they need to share 
information using wireless communication. 

Even though the dynamic model of the UAVs when scan¬ 
ning a target can be highly nonlinear and possibly complex. 


where the state vector Xj(f) is a random vector in K", which 
can contain the velocity, acceleration or angular velocity, etc., 
and it depends on the way the mobilities of nodes are mod¬ 
elled. The additive uncertainties w.;(f) £ R n (assumed to be 
a second-order moment process) can represent the airflow for 
aircrafts (see Section |IIJ, velocity uncertainty for mobile phone 
users, etc. The location of node i in the d-dimensional space is 
denoted as y i(t) £ The constant matrices A, £ M nx ", and 
C i £ M rfx " are the model parameters determined by the node 
dynamics. The initial vector for the differential equation |2]) 

2 Simulation Parameters: Signal power is 1 for each UAV, and the mobility 
model follows Q with too = 0 J 2 = O.lrad/s and = — O.lrad/s. The 
initial locations for UAV 0 - 2 are (500, 500), (-400, -300), (400, 0) 
(selecting the target center as the origin), thus according to jT| average radii 
are Rq = 500\/2m, Ii\ = 500m, and R% = 400m. The powers of ] ^ (t) 
and w[ 2 \t) are assumed to be unit. We assume the Path loss function r ~ 2 , 
where r is the Euclid distance between transmitter (UAV) and the point whose 
interference is needed to be calculated. 
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is Xj(s) at time s and is independent from vrfit), fit > s, 
because wfit) only affects the future behavior of Xj(£). 

For the UCM example in Section [II] we have 


A* = 


' 0 


0 > 


C, 


ri oi 

o 1 > 


(4) 


and indeed yfit) = x, (f). This shows UCM is a second-order 
linear mobility model. 


Definition 1 (Homogenous Mobility). If pairs (A,. C,) for 
all the nodes are equal, i.e.. 


Aj = A, C,=C, *G IV} (5) 


and wfit) /or all nodes are i.i.d., then the node mobilities are 
homogenous. 


Remark 1. The random walk and discrete-time Brownian 
motion models given in are homogenous and can be 

regarded as a special case when 0 is discretized and A = 0. 
Thus, they are homogenous first-order linear mobility models. 
The one-dimensional homogenous continuous-time mobility 
model with 


ro l 

0 ln(l — a) 


C = [1 0], w (t) 


0 

u^ 2 fit) 


( 6 ) 


can be discretized to recover the Gauss-Markov mobility model 
given in [20], where a € (0,1) and the mean of u^ 2 \t) is 
the asymptotic velocity. 


Remark 2. For any second-order moment process, the mean 
vector E[y.;(t)] and covariance matrix Cov[yj(£)] of y fit) 
can be calculated, even though the pdf of y fit) has no 
closed-form expression in general. For a fixed covariance 
matrix Cov [y,: (£)], the entropy is maximized when y fit) 
is Gaussian © It implies that for all kinds of second- 
order moment process w fit), Gaussian process contains the 
most uncertainties. Practically, it is difficult to determine 
what kind of process w fit) is, and hence the best choice 
is to conservatively (since it contains the most uncertainties) 
regard w,; (t) as a Gaussian process. Furthermore, w fit) is 
independent on time t in practice, we often consider it as 
a white noise. Therefore, wfit) can be regarded as Gaussian 
White Noise ( GWN). Actually, GWN is widely used in modeling 
uncertainties like disturbances I\2lfi - f28^j . In the rest of this 
paper, we investigate the statistics of y fit) when ~wfit) is 
GWN. It can be proved that, if w fit) is Gaussian in <0- then 
y fit) in 0 is still Gaussian. Thus, the pdf of y. t (t) can be 
determined by its mean and variance. 

Lemma 1 (PDF of Node Distribution). Assume that w fit) is 
GWN, the pdf of y fit) at time t is a Gaussian distribution 
with parameters 

E[y,(f)] = Cie Al(t - s) x t (s), (7) 

Cov[ yi (t)] =CiQ Xi (t)Cf, (8) 

where 

®xfit)= [ e A i( t-r ) Cov[w.;(t)] e A ^ t ~ T ' ) dr. (9) 


B. Dynamic Reference Point 

A static reference point can be used when we investigate 
the interference statistics at a fixed point, like at a fixed 
base station. However, if we want to analyze the interference 
statistics of a mobile node in a MANET, e.g., a moving 
robot or a UAV (see Fig. l(b)| >, the reference point should 
be dynamic with possible uncertainties. 

We assume the dynamic reference point, denoted as yo(t), 
satisfies 


x 0 (£) = A 0 x 0 (f) + w 0 (f) (10) 

Yo(t) = C 0 x 0 (t), (11) 

and the relative location of node i from this reference point is 
given by 

yfit) ■= yfit) - yo(£). (12) 

Let f y (yfit)) denote the pdf of y fit). The following Lemma 
derives the pdf of y fit), denoted by f y (y fit)). 

Lemma 2 (Node Distribution w.r.t. a Dynamic Reference 
Point). Assuming the mobility model of nodes and reference 
point are GLC with GWN, the location of the ith node relative 
to the dynamic reference point is Gaussian distributed at time 
t with mean 

E[y/£)] - C i e Ai < i - s >x i ( S ) - C 0 e A ^- s )x 0 ( S ) (13) 

and variance 

Cov[y 4 (f)] = Ci<d x fit)Cj + C 0 Q Xo (t)C^. (14) 

Proof: Lemma [2] follows from Lemma [T] and ( fl2| ). ■ 

IV. Interference Prediction 

The main problem of study in this work is to characterize 
the interference received at a reference point at a future time 
instant given the interferers’ mobility and location information 
at the current time instant, i.e., interference prediction from the 
current time into the future. We use the GLC mobility model 
defined in the previous section to describe the mobility of the 
interferers and the reference node. Specifically, the quantity 
yfit) i e { 0 ,..., N} in Q or © is the random variables 
that describe the locations of nodes (interferers or reference 
node) at time t with a known initial condition at time s, where 
t can be viewed as a future time instant and s can be viewed as 
the current time instant. To make the time-dependency more 
explicit, we rewrite it as yj(£|s) in the remainder of the paper. 

A. Problem Description 

Suppose reference node 0 is receiving information from its 
transmitter. Assume that there are N mobile interferers in this 
network whose mobilities are modeled by ([2] and ([3]), and 
their interference lasts for time duration [to,tf\. The aggregate 
interference on the reference node at time t £ [£ (l , tf], 

conditioned on knowing the interferers’ node dynamics and 
locations at time s with s < t, can be defined as 

N 

i=1 


Proof: See Appendix |A] 


(15) 
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where hi is the multipath fading gain with E[/ij] = 1, which 
is independent of y ): (f|.s) and time t. g(-) denotes the path 
loss function, and || • || is the Euclidean norm. For the general 
case of s < t, the quantity J(f|s) represents the interference 
prediction at a future time instant t based on the information 
available at the current time instant s, and hence, it is a random 
variable due to the uncertainty in the mobility over the time 
duration from s to t. For the very special case of s = t, 
the quantity I(t\t) represents the actual interference received 
at time t which is a constant instead of a random variable. 
In fact, I(t\t) can be viewed as a realization of the random 
variable I(t\s) for s < t. In this paper, we simply refer to 
I(t\s) as interference prediction. 

We consider the problem of predicting the statistics of 
interference received by reference node 0 at time t with the 
available information from current time s < t. Denote 


N 


S[/(t|s)] =S^/i i £f(||y i (f|s)||) 

i= 1 


(16) 


where §[■] can be any statistics of the interference prediction, 
such as the mean value, and yj(f|.s) is a conditioned random 
vector which represents the relative location of interferer i 
from node 0. In most cases, evaluating §[/(f|s)] requires the 
pdf of y l (t\s), which can be determined by Femma[2] 

In the remainder of this section, the mean value and MGF of 
the interference prediction will be considered in Section [IV-B 
and Section |IV-C[ respectively. In Section |IV-D| we will 
propose a compound Gaussian point process functional as a 


general framework to study these statistics. In Section IV-E 
the information decay in interference predictions will be 
discussed, and the prediction has close ties to BPP modelling 
when time goes to infinity. 


B. Interference Prediction Mean 

The mean of the interference prediction is given by the 
following theorem. 

Theorem 1 (The Mean of the Interference Prediction). The 
mean of the interference prediction E[/(f|s)] is 

N r 

E [ / ( i l s )]=H / 9{\\yi(t\s)\\)fy (y l (f|s))d(y i (f|s)), (17) 

i=l J * d 

where f y {y, (t\s)) is the pdf of y i (t\s), and can be obtained 
from Lemma [2] 

Proof: The mean value of interference prediction is 

N 

E[7(f|s)] =^E[J i (f|s)]. (18) 

i=i 

According to ( fl~5j ), E[/i ? ] = 1 and the independence of h., and 
y,j(f|s), we can derive 

E[li(t\s)]= f 5 (||y i (i]s)||)/i,(y i (t|s))d(y i (i|s)). (19) 
JR 11 

Thus, 0 can be obtained. ■ 

The calculation for the mean of the interference prediction 
will be discussed in Section HV-DI 


C. Interference Prediction MGF 

Similar to the mean interference, the MGF of the interfer¬ 
ence prediction can be derived in the following theorem. 

Theorem 2 (The MGF of the Interference Prediction). The 
MGF of the interference prediction is 


E[e^ / ( t l®)] = /E; l [e /3?li9 (ll yi ( t l s )ll)]/ s/ (y i( f|s ) ) d (y i (<|s)), 

i=i J * d 


( 20 ) 


where / y (yj(f|s)) is the pdf of yffs), and can be obtained 
from Lemma [2] 


Proof: With ( p~5| ) and independence of hi and y ft). 
Theorem [2] can be proved. ■ 

If the power fading is Nakagami-jn. (m = 1 gives the 
Rayleigh fading), i.e.. 


fh(x) 


m m x (m—l) e ~mx 

F(m) 


( 21 ) 


then the MGF of the interference prediction can take a more 
specific form stated as follows. 

Corollary 1 (The MGF of the interference prediction 
w.r.t. Nakagami-m Fading). With Nakagami-m Fading m 
E[e /3/ (*l s )] in ( |20| ) becomes 



m 


Vm - j 05 (||yj(t|s) 


fy (y,:(i|s)) d (yi( i |s)), 


( 22 ) 


where /35(||y. i (f)||) < 1. 

Remark 3. Av already discussed, 7(f|s) is a random variable 
that represents the interference prediction at a future time 
instant t based on the information available at the current 
time instant s. The uncertainty of the random variable can be 
computed using its MGF. For instance, one can compute how 
much the realizations of I (t\s) deviates from its mean value 
using the variance 

Var[7(f|s)] = E[I 2 (t\s)} - (E[/(f| s )]) 2 , (23) 

where the first and second moments of I(t\s) are used. Since 
the actual interference received at time t, i.e., I(t\t), is a 
realization of the interference prediction I(t\s), the variance 
computed above tells on average how much the actual inter¬ 
ference received at time t deviates from the predicted value at 
time s using the mean prediction. 

The calculation for the MGF of the interference prediction 
is discussed in Section HV-DI 

Although the expressions of either the mean or MGF of the 
interference prediction do not usually admit any closed form 
due to the generality of the GFC mobility model, exceptions 
are found in some special cases where closed-form expressions 
are obtained. These results are discussed in Remark 0 in the 
next section. 



















6 


D. Compound Gaussian Point Process Functional and its 
Series Form 

In both Section |TV-B| and Section [TV-C| the mean and MGF 
have similar integrals to evaluate. This suggests there may 
be some general methods to compute these quantities. Here, 
we define the Compound Gaussian Point Process Functional 
(CGPPF) as a general framework for computing these statis¬ 
tics. Its series form is also given, which will also be useful 
in analyzing a limiting property of interference prediction in 
Section IIV-EI 

Definition 2 (Compound Gaussian Point Process Functional). 
The CGPPF is a functional G: V-^Mo/ the form 

G[v\=Ey[v(\\y\\)]= f J'(llyll) fy(y) dy, (24) 

JR'* 


where v £ V is a Lebesgue Integrable function, and f y (y) is 
a pdf of d-dimensional Gaussian distribution given by 


f y ( y) 


1 c -4(y-M) T s 1 (y— m) 
(2tt)5|S|5 


(25) 


where p, £ and S £ R dxd are mean vector and covariance 
matrix of location vector y. For example, fi, = E[y,(t|s)] and 
S = Covfy.^tls)] are mean and covariance ofy i (t\s). 

Remark 4. For interferer i and dynamic reference point 0, 
the following are derived: If v{ || • ||) = 5 , (||y i (t|s)||), then (|24]) 
returns the mean of the interference prediction from interferer 
i., Ifu(\\-\\) = E/j [e^(t)9(lly,(fis)ll)] ; then <[24]» gives the MGF 
of the interference prediction from interferer i. Note that the 
path loss function, g(-), can be arbitrary. 

Remark 5. In most cases, this functional cannot be simplified 
to a closed-form expression, and numerical integration is 
needed. Nevertheless, if 


p = 0, S = diag{cr,..., cr} (26) 


is satisfied, we can get closed-form expressions for the first 
and second moments of the interference prediction, where 
cr > 0 is the std (standard deviation) of all components in 
location vector y. These expressions are given in Appendix] D\ 
The usefulness of condition ( |26| ) will be further discussed in 
Section \IV-E\ 


If the integral in ( |24| ) exists, the CGPPF can be expanded 
into a series form, which is the cornerstone for analyzing the 


limit properties of interference predictions in Section IV-E 


Theorem 3 (Series Form for Compound Gaussian Point 
Process Functional). Let P be an orthogonal matrix such that 
P T P = V = Wi}dx l and P T £ _1 P = diag{l/cr},..., 1/crJ}, 
then G[v] can be written as 


1 

(2tt)5|E|5 


lim V 

R—foo * ^ 
n—0 


(-i)" 

2 n n\ 


E 

ki+k 2 +k 3 =n 


n 

fci, fc 2 , h 




( 27 ) 


The parameters of © are listed as follows: 
r R 


'PM = 


■(r) r 2fc i +A:2+d - 1 dr, 


/o 


0 =[£(*)' 

a= 1 


2-1 k 3 


E 


4 i ) + ... + 4<*) = 

4 1) +-+4 d) = 


fci 

r.(!) i.( d 

"a , - - -, 


z.C 1 ) u 

"'2 >■■■■> ^2 


(d) 


(28) 


„ (29) 


n f- 


1 <a<d 
1 <b<d 


(30) 

Jj^ a ) 2k ^^b) k ^V((f>) d^j, 

where the integration is with respect to the d— 1 dimensional 
vector 4> = (f>i, 4 > 2 t ■ ■ ■, 4>d-i over the domain © = {</>: 
4> £ [ 0, 7r] x • • • x [0, 7r] x [0, 27t]}, and 


d—2 


$ d ] T = 


cos fi , sin fi cos f> 2 , • ■ •, 

d -2 d—l T (31) 

JJ sin (j) p cos 1 , s in <t> P , 


1 


p =1 


d-1 


V{(j))d(j) = JJ(sin</) Z ) d q 1 d(/) q . 

q =i 

Proof: See Appendix [B| 


(32) 


Remark 6. It should be noted that there are two required 
integrations in Theorem [d] For ( |28[ ), it has a closed-form 
expression when calculating the mean and MGF of the in¬ 
terference predictions with the widely-used path-loss function 
of the form 


s(llyll) 


i 


(33) 


where e > 0 (here e = 0 refers to a singular path loss), and 
a is the path-loss exponent. The expression of ( |28| ) is given 
in Appendix [c] 

To evaluate in ( |29| ), we need to compute the integral 
in P0[ ), i.e., 

f ($ a ) 2fc ^($ b ) fe ^V(0)d0, (34) 

J® 


the calculation of which is rather complex. However, in the 
real world, 1 < d < 3, and it is relatively easy to deal with. 
In Appendix\C\ we give the closed-form expression of ( |29| > for 
the cases of d = 1 (e.g., a vehicular network on a highway) 
and d = 2 (which is the most common scenario of interest). 


E. Gaussian BPP Approximations for Interference Prediction 

The interference prediction method discussed in previous 
subsections is naturally based on the mobility model of 
individual interferers. Although prior work on interference 
analysis for MANETs did not explicitly study the interference 
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prediction problem, one useful idea from them is to use a 
certain point process to approximate the node distribution 
in the network from which the time-invariant interference 
statistics can be derived. For networks where its node dis¬ 
tribution (in distant future) can be well approximated using a 
point process, the information on the initial positions of nodes 
and node mobilities becomes unnecessary in determining the 
interference statistics. Clearly, not all mobile networks can 
have such a nice property. In this subsection, we study the 
condition under which the interference prediction at a time of 
far future can be well approximated based on a simple point 
process. We will see that the series form of the CGPPF will 
be useful in determining such a condition. 

Firstly we give the definition of Gaussian Binomial Point 
Process (Gaussian BPP). 

Definition 3 (Gaussian BPP). Let f y be a pdf of Gaussian 
distribution with support R d . A Gaussian BPP with N points 
on R d is a set of i.i.d. random vectors {yi,..., yjv}> each 
with pdf f y . 

Due to the effect of Gaussian white noise, i.e., w fit) in |2|, 
the uncertainty of our prediction will increase with time, which 
implies the information available for prediction will decay. In 
this section, we will focus on the prediction into the far future 
with N interferes governed by homogenous mobilities. The 
following theorem gives a necessary and sufficient condition 
that the statistics of interference predictions can be approxi¬ 
mated as those generated by interferes whose locations follow 
a Gaussian BPP when t becomes large enough. 

Theorem 4 (Necessary and Sufficient Condition for Gaussian 
BPP Approximation). Assume that all interferes’ mobilities 
are homogenous, as defined in Definition [7] and the integral 
in exists, Vi/ € V, the predictions satisfy 


lim (Gfiy\ - Gj[v\) = 0, Vi, j £ ,N} 


t—f OO 


(35) 


if and only if 


Jf) JJ) 

V 'la r 'la 

lim —TT7 — lim —rrr 

t—¥ oo _vv t—foo 

On O a 


= 0, 


Va = {!,.. .,d), Vi, j £ ,N} (36) 


lim [Gfiv] - Gfiv}} 


can be written as 


1 

(2tt)5|S 1 |5 


lim > 

R—f oo ^ J 
n—1 


(~i) ra 

2™n! 


E 

fcl+fe2+ fc 3 = ri 

fc 3 ^0 


n \ 
ki,k 2 , k 3 ) 


lim (fl; 

t-S-oo V 




(38) 


Obviously, formula ( [38] ) equals 0 if and only if ( [35] ) holds. 

Necessity. By the contrapositive, if the ([36]) does not 
hold, lim (f2j — ilfi fi 0 (since arbitrary v determines 
arbitrary ( |34[ i in ( |3()[ . if ^lim ifl, f> 7 ) = 0, ( |36[ i must be 
satisfied). As a result, the ( |38[ > does not equal 0. Therefore, 
the contrapositive is established, and the necessity of ( [36] ) is 
proved. 

Sufficiency. If ( [36]) is satisfied, then the ( [38] ) equals 0, so as 
for (|35|). Thus the sufficiency of (|36|) is established. ■ 


Remark 7. Theorem [7] tells us the interference predictions, 
from interferes with homogenous mobility asymptotically con¬ 
verge to the predictions from a Gaussian BPP. Nevertheless, 
i/m does not hold, we cannot use the BPP approximation. 
In Section [V] we will see examples in both scenarios. It should 
be noted that rja ^ and ain ( |36| ) can be easily calculated 
for a given GLC mobility model (see Section [F[/or examples). 


Remark 8. Gaussian BPP approximation implies that the 
mean and MGF of the interference prediction become 


E[/(f| S )] 

E[e /3/(t|s) ] 


N G l [s(||yi(f|s 
G 


Vi 


IAm-/3s , (||y i (f|s)[|) 


N 


Vi, 


(39) 

(40) 


when t y>> s. Therefore, we can calculate the mean or 
MGF of the interference prediction using just one (arbitrary) 
interferer’s information instead of all interferes, and the 
computation is significantly simplified. Note that ( |39| ) and (gOj) 
change with time t, because the distribution of the predicted 
location y(f|s) varies with t. 


Corollary 2 (Necessary and Sufficient Condition for Gaussian 
BPP Approximation with pt = 0). Assume that all interferes’ 
mobilities are homogenous and the integral in m exists, Vi/ £ 
V, the predictions satisfy 


lim (Gi[v\ - G 0 [f) = 0, Vi £ {1,..., N} 


(41) 


holds, where rja 1 and aare defined in Theorem [i] with 
superscripts (i) for ith interferer and 


Gi 


|y»ll) 


i 


(27r)f|Si|3 


s -|(y 1 (y i -fi i ) ( jy i _ 


(37) 


In ( |37| ), pti, S i and y i stand for E[yfit\s)], Cov[yj(f|s)] and 
yft\s), respectively. 


if and only if 

(*) 

, lim ~Ti) = Va £ (1,..., d}, Vi £ {1,..., N} (42) 

t—fOO rr \ 1 ') 
v a 

holds, and 



where pt = fi 1 = ... = pL N (fi i = IE [y^ (t| s)] ^ and S = 
Ei = ... = Sat (Sj = Cov[yj(f|s)]). 


Proof: If the mobilities of all interferers are homogenous, 
from Lemma [2] the covariance matrix E, in <HZ} for all the 
interferers are the same. Thus, if the integral in (|24|) exists, 


Remark 9. Corollary [2] implies that we can use the Gaussian 
BPP with p = 0 to approximate the mean or MGF of the 
interference prediction when condition © holds. It should 











be noted that it is easy to test condition © by calculating 
rj'a > and for a given GLC mobility model. Furthermore, 
if X = diagjcr,..., a} also holds, we can use the result 
stated in Remark [5] to obtain closed-form expressions for 
first and second moments of interference prediction given in 
Appendix [S] An example of mobility model that has such nice 
properties is Brownian motion. 


V. Simulation Examples 

In order to corroborate our theoretical results, simulations 
are presented to illustrate the effectiveness of interference 
prediction. We consider three different mobile networks with 
the mobility model given by: (1) the basic 2-dimensional 
Brownian motion, (2) 2-dimensional Brownian motion with 
inertia in the velocity and (3) UCM in 3-dimensional space. 
The basic Brownian motion is a simple and commonly-used 
mobility model while Brownian motion with inertia is a new 
mobility model that has not been considered in the literature. 
We will study and compare the interference prediction results 
for these two examples. Specifically, we will see that the 
existence of inertia completely changes the limiting behavior 
of the interference prediction. After that, we move from 2- 
dimensional examples to a 3-dimensional example of UCM 
which can represent the scenario of UAV target scanning 
similar to that in Section II. Note that all three mobility models 
are special cases of the GLC mobility model developed in this 
work. 

To make predictions, the parameters for path loss function 
(e.g., the path-loss exponent a) and multipath fading (e.g., the 
value of to for the Nakagami-m fading channel) are necessary. 
In addition the following parameters are assumed to be known 
at time s: 

• For 2D Brownian motion (Section | V-A| >, the location for 
each node. 

• For 2D Brownian motion with inertia (Section [V-B| >, the 
location and velocity for each node. 

• For UCM (Section [V-C| i, the location and angular velocity 
for each node. 

Note that the node location and velocity at only one time 
instant (not frequently updated) can be obtained in many 
practical scenarios [29J-|31|. For mobile users (Section V-A| >, 
their locations can be updated by Global Positioning System 
(GPS) (29). For vehicles (Section |WB) or UAVs (Section|WC), 
their locations or velocities can be obtained by Differential 
GPS (DGPS), e.g., (30), (U], and the angular velocity can be 
calculated by the node location, node velocity and the center 
location of UCM. 


A. 2D Brownian Motion 


We consider a network having 7 nodes with mobility 
governed by Brownian motion, which can be used to describe 
human mobility (32). As a special case of our GLC mobility 
model, 2D Brownian motion has the parameters 


Aj 


ro oi 

0 0 ’ 


C* 


n oi 
o 1 > 


(44) 


where i = 0,1,...,6 (that is, 6 interferers and 1 reference 
point). Note that the reference node also does Brownian 



(a) 



(b) 

Fig. 2. (a) Node motions, (b) Aggregate interference at reference node 

and the corresponding predictions. E[/(f|s = 0)] stands for the mean 
value of interference prediction based on information available at s = 0, 
and std[/(t|s = 0)] is the standard deviation of prediction. They are both 
dependent on time t and s. Each error bar E[/(t|s = 0)] d= std [7(£|s = 0)] 
is symmetric about E[/(t|s = 0)], and the negative part of E[/(t|s = 
0)] — std [7(£|s = 0)] is omitted. 


motion. Furthermore, all the nodes (including the reference 
node) start moving from the origin at time to = 0, i.e., 
yz(0) = Xz(0) = 0. 

The uncertainty w, follows 

w»(t) = [w ( i 1 \t),w{ 2 \t)] T , (45) 


where all wfit) for all i are identical GWNs with unit power. 
A realization of Brownian motion is shown in Fig. |2(a)| 
Considering e = 1 and a = 4 in the path loss function and 
Nakagami fading with m = 2, we can predict the mean and 
standard deviation of the aggregate interference at node 0 from 


s = 0, i.e., E[/(f|s = 0)] and std[/(i|s = 0)], see Fig. 2(b) 


The std[/(t|s = 0)] measures how much the realization of 
/(i|s = 0), i.e., I(t\t), deviates from its predicted mean value, 
and thus std[/(f|s = 0)] reveals the uncertainty of interference 
prediction. In this case, both E[/(t|s = 0)] and std[/(f|s = 0)] 
have closed-form expressions (because (26) in Remark [5] is 
satisfied), which can be derived by calculating E[/j(£|s = 0)] 
and std[/i(f|s = 0)], Vi, from © and (78) in Appendix [d] 
From Fig. |2(b)| we can see that the mean of the interference 
prediction does not perform well, and the Coefficient of Vari¬ 
ation (CV) std[/i(f|s = 0)]/E[U(f|s = 0)], i.e., normalized 
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Fig. 3. Predictions on the mean of aggregate interference based on frequently 
updated location information at time s = 0, 4.6, 9.7,14.8,19.9, 25. 


standard deviation, increases with t. This is due to the fact 
that the uncertainties Wj(t) in velocities x t (f) dominate the 
behavior of the mobilities of nodes. In Fig. |2(a)| the node 
locations are far away from the origin, while prediction tells 


us E[y.j(f|s = 0)] = yj(0) = 0 (see Fig. |2(a)[ ). 

In Fig. |2(b)| s = 0 means the information on node locations 
is updated at time 0. If we want a better prediction for near 
future, the information on node locations need to be updated 
more frequently. Interference prediction is then done with 
updated location information until the next update (see Fig. [3]>. 

We end this subsection with a discussion on the Gaussian 
BPP approximations for these predictions in Fig. [3] Actually, 
these predictions can be approximated by a Gaussian BPP at 
origin. This is because the condition ( |42| in Corollary [2] is 
satisfied. 

(i) (i) 

It should be noted that r/„ ' and a a can be easily calculated 
to test the condition ( |42| : By Lemma[2] the mean and variance 
of the location prediction at t, from s for zth node are 


M = Yi(s) ~ yo(s), E = diag{2(f-s),2(f-s)}, (46) 


where y*(s) — yo(s) is finite. As E is naturally a diagonal ma¬ 
trix (which means the orthogonal transformation P T S 'P = 
diag{l/erf,..., 1/cr^} is not required), we have 


Vi = Mi = y ( i 1 \s)-yo 1 \s), r) 2 = y-2 = y^\s)-y^ 2 \s), 

cri = cr 2 = \/2(t - s). (47) 

Thus, the condition ( |36| ) in Theorem [4] is satisfied, which 
implies the predictions made from different s in Fig. [3] can 
be approximated by the corresponding prediction from the 
Gaussian BPP with 6 nodes whose location pdf has /x = 0 
when t — s is large enough. The approximation has the 
following close-form expression 


E[/(f — s|s)] = 6E [Ii(t — s|s)] 

^ 6 Ci (#) sin 2 ^+ 3cos # 

2y / ecr 2 

where E[/,(£ — s|s)] is derived from ( |77[ ) in Appendix [ d] and 
a = ai = (72 = \J2(t — s). Since E[/,(£ — s|s)] is time- 
invariant, it is exactly the E[J(f|s = 0)]. Therefore, if we 



translate the time-axis by s — t, these prediction E[/(f|s)] 
(s = 4.6,9.7,14.8,19.9,25) will asymptotically converge to 
E[/(f|s = 0)] as t — s increases (see Table 0>- 


B. 2D Brownian Motion with Inertia 


The basic Brownian motion in Section V-A is dominated by 


velocity uncertainty, thus our predictions cannot be expected 
to offer great accuracy. In this section, we focus on the 
Brownian motion with velocity inertia in 2D space, and it can 
be employed to describe the motion for vehicle or pedestrian 
with destination. 

The mobility model has the parameters 


(49) 



ro 

1 

0 

01 






A, = 

0 

0 

0 

0 

0 

0 

0 

1 

, Q = 

ri 

0 

0 

0 

0 

1 

°p 


0 

0 

0 

0 







where i = 0,1, 


where x^(0) = 0 and x ( f J> (0) = 0. Note that the zth node 
location is given by yz(0) = [x- 1 ^(0),a;- 3 ^(0)] T . The non¬ 
zero states, a^(0) and ^^(O), which represent the zth node 
velocity components (i.e., yz(0) = [a^ 2 ’ ) (0),a;^ 4 ' ) (0)] T ), are 
randomly generated in [—1,1] and assumed to be known. The 
W ? ; follows 


, 6, and with the initial condition 


( 0 ) = [ 0 , 4 2 ) ( 0 ), 0 ,^ ( 4 ) ( 0 )] 


(50) 


„(3) 


W. 


(t) = [wf } (£),0] 


(51) 


where all 1 (t) are identical GWN with unit power, which 
implies the node velocities are fluctuating due to these uncer¬ 
tainties. The parameters of the path loss function and multipath 


fading are the same as those in Section V-A (i.e., e = 1, 
a = 4, zrz = 2). A realization for Brownian motion with 
inertia is shown in Fig. |4(a) The mean and standard deviation 
of interference prediction from s = 1.8 are shown in Fig. |4(b)[ 
From Fig. 4(b) we can see that the prediction performs much 
better than that for Brownian motion without inertia, and the 
standard deviation is smaller than that in Brownian motion 
without inertia for the same interference level. It also should 
be noted that, similar to the Brownian motion without inertia, 
the CV becomes larger with increasing t. 

In terms of the limiting behavior of the prediction, unfor¬ 
tunately, it cannot be approximated by the prediction from a 
Gaussian BPP in Theorem [4] no matter how large t is. It can 
be calculated that 

m = x { P (s) - 4 1} ( s ) + [ x f ] 0) - 4 2) (s)] ■ (t ~ s) 

V 2 = x^is) - 4 3) 0) + [^ 4) ( s ) - 4 4) (s)] • - s) (52) 

0-1 = CT 2 = \j 2 [t - s). 


Because \_x\ (s) — (s)] ^ \x\ (s) — Xg (s)] at s = 1.8, 

condition ( |36| in Theorem [4] is not satisfied. The prediction 
based on the actual mobility model and that based on BPP 


approximation are shown in Table III and it turns out that 
there is a big difference between them. 
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TABLE II 

Convergence Behavior for the Mean of the Interference Prediction in Networks with Brownian Motion Mobility models 




t — s = 10 

t - s = 50 

t — s = 100 

t - s = 500 

E[I(t\s 

= 0)] 

0.2201 

0.0463 

0.0233 

0.0047 

E[/(f|s 

= 4-6)] 

0.2115 

0.0459 

0.0232 

0.0047 

E[/(fjs 

= 9.7)] 

0.2032 

0.0455 

0.0231 

0.0047 

E[I(t\s 

= 14-8)] 

0.1772 

0.0442 

0.0228 

0.0047 

E[/(fjs 

= 19.9)] 

0.1768 

0.0441 

0.0228 

0.0047 

E[/(fjs 

= 25)] 

0.1850 

0.0446 

0.0229 

0.0047 


TABLE III 

Comparison of the Mean of the Interference Prediction in Networks with BPP Approximation 


t — s = 10 t — s = 50 t — s = 100 t — s = 500 


E[I(f\s = 1.8)] 2.679 x 10' 2 1.837 x 10 -4 3.025 x 10 -5 8.665 x 10 -7 

BPP Approximation 0.2201 0.0463 0.0233 0.0047 




(b) 

Fig. 4. (a) Node motions, (b) Aggregate interference at reference node and 
the mean of the interference prediction. 


1 and 2), whose parameters of mobility models are 


A, = 



' 0 

-0.1 o' 



T 

0 

o' 

A 2 — 

0.1 

0 0 

, 

-0 

= c. 


0 

1 

0 


0 

0 0 



0 

0 

1 

0 

0.1 

o' 


T 

0 o' 




-0.1 

0 

0 

. C 

1 — 

0 

1 0 




0 

0 

0 


0 

0 1 





(53) 


which implies the hover angular velocities for UAVs 0 and 
2 are both — O.lrad/s, and for node 1 is O.lrad/s. The initial 
locations yi(0) = Xi(0), i = 0,1,2 are 


yo(0) = 


'500' 


'-400' 


'400' 

500 

> yi(o) = 

-300 

, yo(0) = 

0 

500 

700 

300 


(54) 


where yg 3 \o) = 500, y^(0) = 700 and y^\o) = 300 are 
the initial altitudes for UAVs 0, 1, and 2. The w; follows 


W,; 


where all wf ' 1 (t), j = 1,2,3 are identical GWN with power 
tj 2 = 100. The parameters of path loss function are the same 
as those in Sections [V-A| and [V-B| (i.e., e = 1, a = 2). Because 
there are few obstacles in airspace, we assume that there is no 
mult ipath fading. The UAVs’ motion-trajectories are shown in 
Fig. 5(a) The aggregate interference and E[/(f|s = 17)] are 
shown in Fig. |5(b)| We see that the interference prediction 
is very close to the actual interference. It is interesting to 
see that the CV is not always increasing with t. However, 
if we consider the interference at the same level, the CV does 
increase with time, e.g., t = 30 and t = 60. 


it) = 1) (t), ^ 2) (t), (*)] 


,( 3 ) / 


(55) 


C. Uniform Circular Motion in 3D Space 

In this section, we revisit the UAV target scanning problem 
discussed in Section [II] In the real world, UAVs that execute 
the scanning task seldom hover in the same 2D plane (their 
flight altitudes differ) for collision avoidance. Thus, UCM 
should be modeled in 3D rather than 2D. Assume there is one 
receiver (reference node, UAV 0) and two interferers (UAVs 


VI. Conclusion 

In this paper, the interference prediction problem in 
MANETs has been investigated. The GLC mobility model 
has been proposed to describe or approximate a large class 
of mobilities in the real world. The statistics of interference 
prediction with respect to a dynamic reference point have been 
analyzed, including mean value and the MGF. We have defined 
the CGPPF as a general framework to compute the mean and 
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- reference node 



interference prediction converges to that from a PPP. Another 
interesting direction for future research is to generalize the 
assumption for w, (t) beyond Gaussian. Last but not least, the 
results on interference prediction obtained in this work can 
also been used to predict the outage probability at a future 
time instant. 


Appendix A 

PDF of Node Distribution 


A more general proof can be found in [261. For complete¬ 
ness, we provide a proof of Lemma [I] here: 

The solution of (2]i is 



Xj(f)=e Ai(t to) Xj(f 0 )+ [ e Ai(t r) Wj( t) dr (56) 

Jt 0 

and then the mean of x,(t) can be derived as 

E[xj(f)] = e Ai ^ _t °^Xj(f 0 ) + [ e A ‘ (t ~ T) E[w;(t)] dr 

Jto (57) 

( = e A *(‘- to )x,(f 0 ), 

where (a) is established by E [w,;(t)] = 0. Then we can get ([TJ) 
from (3). 

The covariance of x, (f) is given by 


(b) 

Fig. 5. (a) UAV locations versus time t, starting from the triangular marks 
and ending at the circular marks, (b) Aggregate interference at reference node 
and the mean of the interference prediction. 


MGF, and discussed important special cases where closed- 
form expressions can be obtained. With expressing the CGPPF 
in series form, we have analyzed the limiting behavior of the 
statistics of the interference prediction, and given the necessary 
and sufficient for when the node locations can be regarded as 
a Gaussian BPP when analyzing the statistics of interference 
prediction. 

The presented work serves as the first step to develop 
more comprehensive results in interference predictions. Even 
though the CGPPF provides a general framework for calcu¬ 
lating the statistics of interference prediction, the closed-form 
expressions exist only in some special cases. It is possible to 
derive closed-form approximations for these statistics using a 
cumulant-based approach (e.g., similar to (33), (34)). Further¬ 
more, for the limiting behavior of interference prediction with 
homogenous mobilities, it would be desirable to obtain a more 
direct link between the mobility model and the condition for 
the BPP approximations stated in Theorem [4] For example, 
we observed from our numerical results that if the matrix A 
of the mobility model is Lyapunov-stable (24) , the condition 
for Gaussian BPP approximation holds, and vise versa. On 
the other hand, it would be interesting to extend the results 
on finite number of nodes to infinite number of nodes and 
study the condition under which the limiting behavior of the 


Cov[xj(f)] ( = Cov[e Ai(f to) Xi(f 0 )] 


+ Cov 


e A 4 t T ^w i (r)dr 


Jto 


(58) 


= 0 + Cov 


e A; ( t r ' l Wi(r)dT 


where (a) follows from the independence of Xj(io) and 
Wj(r) (r € [fo,f])- Additionally, 


Cov 


pA 


e~~"~ w Wi(r)dT 


'to 


= Cov 

rt rt 


e A i(i-r) 


r *0) dr, 


C A iit-ip) 


'to 


’ r 'w i (V’)dV’) 


f [ Cov(e 

J to Jto 


A i(t-r) 


'to Jto 
rt rt 


/j(r), e A d t -V') w ;W’)) dr dtp 


a A i(t-T) 


Cov(w,;(t), Wi(if>))e Ai (i 


= e Xi (f), 


where (a) follows the properties of GWN. Then (8]i can be 
derived from (3]). 


Appendix B 

Proof of CGPPF Series Form 

With the orthogonal transform z = P T y such that P 7 /', = 

V = [Va]dx l and P 7 ’S _1 P = diag{l/cr?,..., 1/crJ} = A, 
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equation \2A\ can be rewritten 


as 


GW] = / Kllyll) 


1 


R d (27r)2|S|2 

1 

dT, J 


0 -s(y-ri) T s 1 (y— m) 


dy 


= / V z 


/ R<* (27r)2|S|2 

( a) i r 


,-l( z ” ? t) rA ( z ~ ? 7) d z 


(2-7r) 2 |S| 2 Jo 


(r)- 


e -i(r*-r7) T A(r#-?7) r ,d- 1 


dr / V(4>)d4>, 


7(r,$) = ^ 


n—0 


2 n n\ 


|z2(r)i 


d-1 


[(r<j> _ rj) T A (rd? — ry)]” |dr. 
By employing the multinomial theorem, 


[(r$ - rj) T A(r$ - T ]))" = 

&l+^2+' 

■51 


n 

ki,k 2 ,k 3 


k 1 +k 2 +k 3 =n . 


where 


E 




Va=l . 


E 


k[ i) +-+k{ d) =k 1 


d / -r \2k± 

n,e 


ki 

hi 1 ) ui d 

Aii , . . . , rv-^ 


E- 

^&=i 


\k: 

2\ 


E 


n 

6=1 


-|-/c^=fc2 

2 


kC 1 ) L,( d ) 

™2 5 • • • ? ^2 


2^6 Y 


t- 2 / 


Thus, ( |59| ) can be written as {27}. 




rt c 


(t+c)e 


c-t-a—1 


where c = 2fci + fc 2 + d— 1 and H 2 F\f) is the hypergeometric 
function (35). When e = 0, the condition c — a +1 > 0 should 
be satisfied, this is because the singularity at 0. 

In order to calculate the MGF of interference prediction, we 


(59) 


where (a) follows a d-dimensional spherical transform. $ = 
[3vUxi and V(4>)d<f> are shown in ( |3 1 [ i and ( |32) . respectively. 
Labeling 

poo 

V(r,&)= / (60) 

Jo 

we expand it into the Taylor Series 

(—l) n r°° 


set v (r) = 

^W] = RC+1 ( e-/3 + R a r m 

tt n (i l+c + a(l-m) 1 + c + a 

R 2*1 1,-j- ; 

\ a a 


thus ^[ 12 ] is 

e (e — /? + R a ) 

(1 +c)(e-/3) 

R a 


P - e 


R l 


( „ \ 7 


H 2 F 1 


( 1+c+q l+c+2q it Q 




/5-e , 


1 + c + a 


( 66 ) 


where e ^ j3. 

B. Derivation for i l in 
For d = 1, 

O = 


)_2) fc 2^fc2 + 2fe 3 


(61) 


(67) 


(j‘2(k-\ J rk 2 J rh 3 ) 

In order to simplified the discussions for the cas d = 2, we 
define two functions 


EEo r/ 2 )(-i 


0 

v^r(^±i) 


(62) 


r (- 


m is even, m > 0, n > 0 
m is odd, m > 0, n > 0 
m = 0, n > 0 

( 68 ) 


and 


I[o,2n](fn,n) = 

T(to, n) n is even, m > 0, n > 0 
0 n is odd, m > 0, n > 0 ; (69) 


) m > 0, n = 0 


(63) 


where 


1=0 


Therefore, with d = 2, we have 

k 3 


(64) 


Q = 


where 


'2 , \ 2 

S’ 

. 2=1 


E 

k^+kif^=ki 
fci 1 - 1 +fc( 2 ' =k 2 


ki 

1.(1) 7.(2) 

Ttl , h, 1 


k 2 

k {1) k (2) 


(70) 


Appendix C 

Integrations in CGPPF Series Form 
A. Derivation for [v\ in {28} 

In order to calculate the interference prediction mean, we 
set v (r) = e+ 1 r „ , thus \D[i/| is 

t R c+1 H 2 F 1 (l.2+l. c + 1+q Hf.) 

-1-2-“- e_J_ £ > Q 


n 

l<i<2 

1<J<2 


2k 


(0 


uO) 


7 2 l. 


7[o,27t] (2/c{ 1} + ^ 1} ,o) j [0)2w] (2fc«,/cf 


e = 0,c — a + l>0, 
(65) 


fas*] J [0i2ot] (0,2^ 2) +fc, 


(71) 
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Appendix D 

Some Closed Forms for CGPPF 


If condition ( |26] > is satisfied, some closed-form expressions 
for E[7j(f|s)] and E[/ 2 (i|s)] can be derived. Please refer 
to ( |33j ) for the definition of a and e. 

Firstly, if e = 0, which implies a singular path loss, we can 
derive the closed-form expressions for first and second order 
statistics as 


E[/i(i|s)] 

E[i?(i|s)] 


r v d r d ~ i v d r[^\ 

/ - 2 - e 2 * z dr= —-— d 

Jo ( 2 tt )2 a d r a cr a 

m +1 f°° Vdr^ 1 _ r 2 

- / - d - e ^ dr 

171 Jo ( 27 t )2 (r d r 2a 

?n2“ +1 7T5(J 2a 


(72) 


(73) 


when gamma function l’(-) has finite value. v d is the volume 
of the d-dimensional ball of radius 1, i.e.. 


2t r d / 2 

Vd ~ T(d/2)' 


(74) 


However, if e > 0, closed-forms are difficult to derive for 
general a > 0. We just give the results for a = 2 and a = 4. 
When a = 2 


E[/j(f|s)] 


v d e^e^T{d/2)TC-^^) 
2(2tt)5|S|5 
(to + l)v d (d — 4)r 
16to(27t) 5 |S|3e 2 cr 4 


{ 


2a 2 e 2 ' 


d 

e 2 



(75) 


2 2 a 



(76) 


For a = 4 and d = 2, 


E[/j(f|s)] = 


2Ci ( A ) sin ^ + cos 


2a 2 


V±_ 


2a 2 




2a 2 


7T - 2Si ( & 


i(#) 


4Vecr 2 


E[/ 2 (f| S )] = 


(to + ljG^g ^0,11/2,3/2 I 16cr 4 ) 


4to£ 2 s / ltCF 2 


(77) 

(78) 


where Si(-) and Ci(-) are sine/cosine integral function with 
the form 


. f z sin a; /'°° 

Si(z) = / -da;, and Ci(z) = — 

Jo x J z 

and G™ n '"’“ p |z^) is the Meijer-G function 
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